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Abstract 

We study the radiative decay within the framework of a phe- 

nomenological approach in which the contributions of cr-meson and w-meson 
intermediate states and of the pion-loops are considered. We conclude that 
the cr-meson amplitude provides the dominant contribution. 
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The radiative decays of neutral vector mesons into a single photon and a pair of neu- 
tral pseudoscalar mesons have been a subject of continuous interest theoretically and are 
becoming an area of active experimental research with the advent of high-luminosity, low 
energy e + e~-machines such as Novosibirsk and Frascati that will make possible studying 
these rare decays with very small branching ratios. Such measurements may serve as tests 
for the theoretical ideas about the interesting mechanisms of these decays because of the 
absence of charged particles hence the bremstrahlung radiation, as well as sheding light on 
the structure of the intermediate states involved in these decays. 

The decay p° — ► 7r°7r°7 among other radiative p-meson decays was first studied by Singer 
IJ assuming that this decay proceeds through an (wir) intermediate state as p° — > (w)+n° — > 
(7r°7) + 7T°. Renard |2| later considered cx-meson intermediate state in addition to this 
intermediate vector meson dominance (VDM) mechanism and noted the dependence of the 
decay rate as well as the photon spectrum on cx-meson mass. 

The radiative decay processes of the type V° — > P°P°>y where P and V belong to the 
lowest multiplets of vector (V) and pseudoscalar (P) mesons were studied by Fajfer and 
Oakes || using a low energy effective Lagrangian approach with gauged Wess-Zumino terms. 
They neglected scalar meson contributions and obtained the branching ratio as B(p° — > 
7r°7r°7)=2.89 x 1CT 5 . The contribution of intermediate vector mesons to the decays V° — > 
P°P°7 was later considered by Bramon et al. [|J using standard Lagrangians obeying the 
SU(3)-symmetry. Their results were in agreement with those of Renard when they 
were corrected using the then present-day data, however they disagree with the numerical 
predictions of Fajfer and Oakes || even if the initial expressions for the Lagrangians were the 
same. Their results for the decay rate and the branching ratio were r(p° — > 7r°7r°7) = 1.6 KeV 
and B(p° — ► 7r 7r°7) = l.l x 10~ 5 , respectively. They also noted that final-state interactions 
could lead to a larger value for the decay rate r(p° — ► 7) through the chain p° — ► 
7r + 7r~7 — > 7r°7r°7. Bramon et al. |J lated considered these decays within the framework of 
chiral effective Lagrangians and using chiral perturbation theory they calculated the decay 
rates for various decays of the type V° — > P°P°7 at the one-loop level, including both 
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7T7T and KK intermediate loops. They showed that the one-loop contributions are finite 
and to this order no counterterms are required. In this approach, the decay p° — > 7r°7r°7 
proceeds mainly through charged-pion loops and Bramon et al. H obtained the contribution 
of charged-pion loops to the decay rate as r(p° — > 7r 7r°7)=1.42 KeV which was of the same 
order of magnitude as the VDM contribution calculated by Bramon et al. || earlier. They 
also noted that the charged-kaon loops gave a contribution to this decay rate which was 10 3 
times smaller than that of charged-pion loops. They obtained the global decay rate and the 
branching ratio by considering the VDM amplitude and the charged-pion loop amplitude as 
r(p° — > 7r°7r°7)=3.88 KeV and B(p° — > 7r°7r°7)=26 x 10~ 6 , respectively, and moreover noted 
that the photon spectrum peaked at high photon energy E T 

Recently, radiative decays of p° and mesons were investigated by Marco et al. || using 
the techniques of chiral unitary theory developed earlier by Oiler et al. j7|. They noted 
that the energies of the two meson systems in the final states of these decays are too large 
for these decays to be treated with standard chiral perturbation theory, and using a chiral 
unitary approach they included the final state interaction of two mesons by summing the 
pion-loops through the Bethe-Salpeter equation. They obtained the branching ratio for 
the decay p° — > 7r°7r°7 as B(p° — > 7r°7r°7)=1.4 x 10~ 5 , and noted that this result could 
be interpreted as p° — > 7<7(7r°7r°) since 7r°7r° interaction is dominated by the cx-pole in 
the relevant energy regime of this decay. Furthermore, they also pointed out that if the 
product of the coupling constants F v and Gy is chosen as negative F v Gy < 0, where the 
coupling constant Gy is deduced from p° — > 7r + 7r~ decay and the coupling constant F v from 
p° — > e + e~, then the result for the branching ratio becomes B(p° — > 7r°7r°7) = 1.0 x 10~ 4 ; 
therefore the measurement of this branching ratio may serve as a test for the sign of the 
FyGy product. 

In this work, we follow a phenomenological approach and attempt to calculate the decay 
rate for the decay p° — > 7r°7r°7 by considering the <r-pole and w-pole diagrams as well as 
the charged-pion loops. We use the coupling constant g po - 7 that we estimated || using the 
experimental value of the branching ratio for the radiative decay p° — > 7r + 7r~7. Since the 
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coupling constant g pal was estimated in terms of the cr-meson parameters M CT and T a) we 
obtain the decay rate r(p° — ► 7r 7r°7) as a function of these parameters and for both signs 

of §p(T7- 

Our calculation is based on the Feynman diagrams shown in Fig. 1. We show the 
amplitude for cr-meson intermediate state in Fig. 1(a). We describe the pcrj- vertex by the 
effective Lagrangian |§ 
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which also defines the coupling constant g prJ1 - For the aim- vertex we use the effective 



Lagrangian [10 
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The decay width of the cr-meson that follows from this effective Lagrangian is given as 
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In our calculation of the invariant amplitude, in the cr-meson propagator in Fig. 1(a) we make 
the replacement M a — > M a — jiT a , where is given by Eq. (3). Since the experimental 
candidate for cr-meson /o(400-1200) has a width (600-1000) MeV |TT|, in our previous work 
||, using the experimental value of the branching ratio for the radiative decay p° — > 7r + 7r~7 
we estimated the coupling constant g pcr7 for a set of values of cr-meson parameters M a and 
T a . In our present work, we use these values of the coupling constant g p(J7 . The upn- vertex 
in Fig. 1(b) is described by the effective Lagrangian |L2] 



(4) 



The coupling constant g^p-n can then be determined using Vector Meson Dominance and 
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current-field identities [|12 



as g t 
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84. The uo'K'y- vertex in Fig. 1(b) is described by the 
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C 



int. 



-Tl-g^^d^daApTT 

3 
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The coupling constant g ttyK1 is then obtained from the experimental partial width [12| of the 
radiative decay uj — > as = 3.315. The p7T7r-vertex in Fig. l(c-e) is described by the 
effective Lagrangian [ |i~0[ | 
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The experimental decay width |ll]] then yields the value = 2.91 for the coupling constant 
g p7rn . We describe the 7r 4 - vertex in Fig. l(c-e) by the effective Lagrangian 
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We note that this effective interaction results in only isospin 1=0 amplitudes. The small 
1=2 amplitudes were also neglected in previous calculations || within the framework of 
chiral Lagrangians utilizing chiral unitary theory. We choose the coupling constant A as 
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and we use — A 
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2 — — as it is obtained in the cr-model with spontaneous symmetry breaking [13 



14. We note that our effective Lagrangians O^Z' and are the ones 
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that result from an extention of the a model to include the isovector p through a Yang-Mills 
local gauge theory based on isospin with the vector meson mass generated through the Higgs 
mechanism |TJJ. Finally, we note that the last of the loop diagrams in Fig. l(c-e) results 
from minimal coupling and ensures gauge invariance of the pion-loop amplitudes. Similar 



loop integrals were evaluated by Lucio and Pestiau |T5[ using dimensional regularization 



and their calculations were confirmed by Close et al. fL6 |. Oiler [17j used their result in his 
calculation of the rate of the radiative decay (f) — > K°K 0, y through an intermediate K + K~- 
loop. In our case the contribution of the pion-loop amplitudes corresponding to diagrams 
in Fig. l(c-e) becomes 
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with a = Mp/M%, b = (p — k) 2 /M%, p(u) and k(e) being the momentum (polarization vector) 
of p-meson and photon, respectively, and 
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In terms of the invariant amplitude .M(E 7 , Ei), the differential decay probability for 



7r 7r°7 for an unpolarized p°-meson at rest is then given as 
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where E 7 and Ei are the photon and pion energies respectively We perform an average over 
the spin states of p°-meson and a sum over the polarization states of the photon. The decay 
width r(p — > 7r°7r°7) is then obtained by integration 
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where now the factor (1) is included because of the 7r°7r° in the final state. The minimum 
photon energy is E 7)m j n . = and the maximum photon energy is given as E liTnax _ = [M 2 p — 
4M^)/2M p =338 MeV. The maximum and minimum values for pion energy Ei are given by 
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We present the results of our calculation in Table 1. Since we use the coupling constant 



g pal that we estimated using the experimental decay rate for the radiative decay p° 
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for a set of values of cr-meson parameters M CT and IV in our previous work || , we also obtain 
the decay rate T(p° — > tv°tv ^) as a function of these parameters. We also consider for a 
meson parameters the values suggested by two recent experiments, that is = 555 MeV 



T a = 540 MeV from CLEO fl§, and M a = 478 MeV T a = 324 MeV from Fermilab E791 
p~9|] . Since in our previous work we obtained a quadric equation for the coupling constant 
g pal which resulted in two values of this coupling constant, one being positive and one being 
negative, for a given set of cr-meson parameters M CT and T a we obtain two values for the 
decay rate p° — > 7r°7r°7, one using the positive value of the coupling constant g pul and 
one using the negative value. We therefore note that the measurement of the decay rate 
r(p° — > 7r°7r°7) may serve as a test for the sign of this coupling constant. This result was 
also noticed by Marco et al. || for the sign of the product of the coupling constants FyGy 
in their work. We further note that the measurement of the decay rate may also help to 
limit the ranges of the cx-meson parameters and IV as well. 

We like to stress that within the framework of our phenomenological approach the main 
contribution to the decay rate T(p° — > 7r 7r°7) comes from the cr-meson intermediate state 
amplitude. Indeed, for Mo-=478 MeV, r CT =324 MeV, and g po - 7 =5.92, we obtain the value 
of the decay rate as T(p° — > 7r°7r°7) = 289 KeV and the contributions of the cx-meson 
intermediate state amplitude, the tc-meson intermediate state amplitude and the pion-loop 
amplitudes are 

r(p° -> 7r°7r°7)a = 323 KeV, T(p° -> vr ^^ = 1-46 KeV, 

T(p° -> 7r°7r°7)^ = 2.7 KeV (13) 

respectively. Our results for the w-meson intermediate state and pion-loop amplitudes are 
in resonable agreement with the calculations of Bramon et al. || and Marco et al. |J. The 
apparent discrepancy between our result for pion-loop contribution and that of Marco et al. 
can be resolved by noting that the coupling constant g pn7T in our work corresponds to 
the constant v f2 p in their work using the framework of chiral Lagrangians. They used the 
values Gy = 55 MeV, f n = 93 MeV resulting in the coupling constant g p7r7r =4.9. On the 
other hand, in our phenomenological approach we determine this coupling constant from 
the experimental decay rate of p° — > 7r + 7r~ [|TTJ using the effective Lagrangian given in Eq. 
6 as g p7T7r =6.0. If we use this coupling constant as g p7r7r =4.9 corresponding to the coupling 



constants in the work of Marco et al. |J , for the contribution of the pion-loops we obtain the 
result r(p° — > 7r°7r°7) 7r = 2.14 KeV which compares well with the result T(p° — > 7r°7r°7) 7r = 
2.11 KeV of Marco et al. [[J. However, since our result for the contribution of the cr-meson 
amplitude is much larger than the contributions of VDM and pion-loop amplitudes, our 
calculation results in a larger value for the decay rate r(p° — > 7r°7r°7) than the previous 
calculations which did not include the contribution coming from the cr-meson amplitude. 

Moreover, if we use the values for cr-meson parameters the values M a =500 MeV, r cr =250 
MeV as used by Soyeur p0| , and furthermore use the value of the coupling constant g pal =2.71 



which Titov et al. || obtained from Friman and Soyeur analysis of p-meson photopro- 
duction using VMD we obtain the decay rate as T(p° — > 7r°7r°7) = 56.2 KeV, or if we use 
M CT =500 MeV, r ff =500 MeV and g prT1 =2.71 we then obtain it as T(p° 7r°7r°7) = 31 KeV. 
In both of these calculations we determine the coupling constant g pn7T from the experimental 
decay rate of p° — > 7r + 7r~ which in the former case results in g p7r7r =4.5 and in the latter case 
in 9(m%— 6-4. We like to stress, however, that in our phenomenological approach, the results 
of which we present in Table 1, we try to estimate the coupling constants we use in our 
calculation utilizing the experimental information, and for the coupling constant g pal we 
again try to relate it to cr-meson parameters through the decay rate r(p° — > 7r + 7r~7) [^]. We 
also use for cr-meson parameters M a and T a the experimental numbers of the fo (400-1200) 
meson ||11|| , which is the candidate for cr-meson, as summarized by Tornqvist pi] . 



The photon spectra for the decay rate of the decay p° — > 7r°7r°7, which may be tested in 
future experiments, are plotted in Fig. 2 as a function of photon energy E 7 . The param- 
eters of the cr-meson adapted for the numerical calculations are the same as above, that is 
M CT =478 MeV, r CT =324 MeV, g pal = 5.92. It is clearly seen that the contributions of the 
cr-term is much larger than the contributions of the w-term, pion-loop term as well as the 
interference terms. The dominant cr-term characterizes the photon spectrum which peaks 
at high photon energies. For comparision, we also plot the photon spectra for the same set 
of cr-meson parameters but for the negative value of the coupling constant g pal = —4.49 in 
Fig. 3. The general shape of the spectrum, its peaking at high photon energy as well as 



the relative contributions of different terms are the same quantitatively as for the case with 
the positive value of the coupling constant g pal shown in Fig. 2. These figures clearly show 
the importance of the interference terms between a-amplitude and pion-loop amplitude and 
also between a-amplitude and ^-amplitude. These interference terms change sign when we 
use negative value of the coupling constant g pal instead of its positive value and this effect 
is also evident in the figures. Therefore, the measurement of the decay rate r(p° — > 7r°7r°7) 
may clarify the nature of the a-meson by limiting the range of its parameters as well as 
determining the sign of the coupling constant g pal which plays a crutial role in analyzing 
photoproduction reactions within the framework of meson-exchange mechanisms such as the 
photoproduction of vector mesons on nucleons near threshold. 
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TABLES 



TABLE I. The decay rate L(p° — > 7r°7r°7) for different cr-meson parameters 



M a (MeV) 


L(j (MeV) 
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L+ (KeV) 
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L_ ( KeV) 


400 


500 


7.66 


5.27 


229 


-5.23 


270 


450 


500 


6.90 


6.19 


256 


-5.33 


243 


500 


300 


4.93 


6.31 


311 


-4.41 


202 


500 


500 


6.36 


7.21 


280 


-5.49 


224 


600 


400 


5.03 


9.28 


334 


-5.63 


198 


555 


540 


6.15 


8.59 


294 


-6.00 


219 


478 


324 


5.29 


5.92 


289 


-4.49 


215 
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Figure Captions: 

Figure 1: Feynman Diagrams for the decay p° — > 7r°7r°7 

Figure 2: The photon spectra for the decay width of p ( 
contributions of different terms are indicated. 

Figure 3: The photon spectra for the decay width of p 1 
contributions of different terms are indicated. 
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Fig. 2 
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